A systematic correlation between two-dimensional flow topology and the abstract 

statistics of turbulence 
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Velocity differences in tfie direct enstrophy cascade of two-dimensional turbulence are correlated 
with the underlying flow topology. The statistics of the transverse and longitudinal velocity differ- 
ences are found to be governed by different structures. The wings of the transverse distribution are 
dominated by strong vortex centers, whereas, the tails of the longitudinal differences are dominated 
by saddles. Viewed in the framework of earlier theoretical work this result suggests that the transfer 
of enstrophy to smaller scales is accomplished in regions of the flow dominated by saddles. 



Two-dimensional (2D) turbulence is a fascinating 
problem with relevance in areas as wide-ranging as the 
dynamics of energy transfer in atmospheric and geophys- 
ical flows [0, ||] to the basic statistical mechanics of inter- 
acting vortices ||, ^, || . Three decades of theoretical and 
numerical work starting from the seminal ideas of Kraich- 
nan |^ and Batchelor Q have provided a picture of 
2D turbulence based upon the scaling laws of abstract 
statistical quantities. This description remains strikingly 
incomplete. For example, there is still no clear physical 
understanding of the mechanisms by which energy and 
enstrophy are transferred between different length scales 
in a turbulent flow; nor is there a conclusive picture of 
how the intense coherent structures that dominate the 
statistics are formed and evolve. 

The challenge is to establish a connection between the 
statistical measures of turbulence and the physical dy- 
namics of the turbulent flow fleld. In the current work, we 
demonstrate that by considering correlations between lo- 
cal flow topology [|[ ^ |l^ and velocity difference 
probability distributions (PDFs) |5[ |l|, 0, || , 
we can make this important connection. The wings of 
the distributions of the transverse and longitudinal ve- 
locity differences are found to be associated with very 
different structures: vortex centers and saddles, respec- 
tively. As a consequence, it will turn out that the transfer 
of enstrophy must be accomplished near saddle points. 
This transfer is the result of a topological asymmetry in 
the turbulent flow manifest in the longitudinal velocity 
difference PDFs. Furthermore, since the wings of the 
longitudinal velocity differences are dominated by strong 
saddles, a complete understanding of intermittency in 
2D must include a motivation for the formation of these 
structures. This powerful technique can be extended to 
other statistical quantities to infer correlations between 
flow topology and turbulent dynamics. 

Local flow topology is characterized by the four first- 
order derivatives in the expansion of the vector velocity. 



The determinant of this Jacobian matrix, 
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A = det(m) = -{lu-^ 
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represents a local balance between the vorticity and 
strain rate. A continuum of structural possibilities exists 
as the relative magnitude of these terms is varied; from 
a symmetric saddle when the strain rate dominates; to a 
linear shear region when the two are of equal magnitude; 
to an axisymmetric vortex center when the enstrophy 
dominates; see Fig. ||a-d. Rivera, Wu, and Yeung 
showed that the probability distribution of the Jacobian 
determinant is non-analytic at A = 0, since centers and 
saddles are topologically distinct, and asymmetric, with 
vortex centers being significantly more likely than saddles 
of comparable strength; see Fig. ^a. 

A statistical correlation between the local flow topol- 
ogy and the velocity differences is constructed from data 
obtained from a 2D flowing soap film experiment, the 
configuration of which is described in Ref. |Q with an 
effective injection scale of 2 cm. The Jacobian determi- 
nant at each location is calculated from the matrix, m, 
averaged over a disk, 57, with center halfway between the 
two velocity measurements and radius, ro, equal to r/2. 
That is, A = det(M), where 
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By performing the average in this manner a scale- 
dependence of the quantity is maintained, allowing the 
method to probe different regions of the enstrophy cas- 
cade. Here, however, we will be reporting results for only 
one separation, r = 0.4 cm. This average over the smaller 
scales in the flow is permissible in the enstrophy cascade 
since enstrophy transfer through a given scale depends 
only on larger scales and not on these smaller structures 

13- . 

There is some sensitivity to flow inhomogeneities in 
this measurement since we are averaging over a macro- 
scopic region of the flow. The inhomogeneity is char- 
acterized by a 25% variation in the turbulence intensity 
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FIG. 1: Fundamental topological structures of 2D incompress- 
ible flow. The Jacobian determinant, A — (1/4)(lj^ — cr^), in 
(a) through (d) reflects a smooth transition from a dominant 
strain rate, a, to a dominate vorticity, uj. 



across the 2.0 cm wide measurement area. This vari- 
ation is due to the particular turbulent forcing mecha- 
nism used (a pair of combs arranged in an inverted wedge 
p2t). Nevertheless, over a 0.28 cm wide stripe down the 
central region of the channel where we gather statistics 
the turbulence intensity varies by less than 1%. Further, 
a simple line average between the two points at which 
the velocity is measured — less affected by cross-stream 
inhomogeneity — produced the same results. 

A rigorous connection between the abstract velocity 
difference statistics and the concrete flow topology is es- 
tablished through multivariate probability distributions 
of the form P{A, Svi(r)), where i represents either the 
transverse, T, or longitudinal component, L, of the ve- 
locity difference. The form of these distributions differs 
significantly between the longitudinal and transverse ve- 
locity increments as a result of the different symmetries 
of the saddle and center; see Figs. |and|. The wings of 
the longitudinal velocity difference PDF are dominated 
by strong saddles, whereas, the wings of the transverse 
velocity difference are dominated by strong centers. To 
understand the reason for this distinction we examine 
the distribution of the velocity differences about these 
two first-order structures in some detail. 

We begin with an examination of a saddle point. The 
symmetry of a saddle is such that the distribution of the 
longitudinal and transverse velocity differences about it 
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FIG. 2: (a) The asymmetric distribution of the Jacobian de- 
terminant, -P(A). (b) Multivariate probability distribution, 
P(A, 5vL{r = 0.4 cm)), indicating the likelihood that a given 
longitudinal velocity difference will be found in a region of 
the flow with topology described by A. Dashed line indicates 
the velocity difference of maximum likelihood. Note that the 
direction of increasing velocity difference along the vertical 
axis is non-standard. Shading and contours represent the log 
of the probability. 



is the same. If the matrix, to, is parametrized as 



X A 



-A 



(4) 



the radial and angular components of the velocity 
about a saddle (described by the symmetric part of to) 
are given by Vr — [A cos(26') -f A sin(20)] r and vg = 
[A sin(20) — A cos(26')] r, respectively. As long as = 
^2_|_^2 j^gy constant, the relative magnitudes of A and 
A serve only to vary the spatial orientation of the saddle. 
We can, therefore, make the simplifying assumption that 
the saddle is described by — A^. Using P{vr{0)) oc 
[deVr{9)Y'^ \e(^^) and P{ve{9)) oc [dgVeiO]]^^ \e{y,) in the 
longitudinal and transverse cases respectively, we obtain 
probability distributions of the form 



P{Sv,ir)\K) 



TT^/r^A^ - Sv^ 



(5) 



regardless of which velocity increment is examined. Here 
we have let, for example, Svi^ — 2vr, and replaced the 
independent variables A and A by an equivalent pair ex- 
pressing the saddle strength, cr^, and orientation. Be- 
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FIG. 3: Multivariate distribution, P{A,5vT{r = 0.4 cm)), 
indicating the likelihood of finding a given transverse velocity 
difference in a region of the flow with topology described by 
A. Dashed lines indicate the velocity difference of maximum 
likelihood. 



cause of this similarity, it's not surprising that the ap- 
pearance of the multivariate distributions, P(A, Svxir)) 
and P{A,SvL{r)), are, at first order, identical for A ^ 
(where A ~ — (1/4)(t^). The velocity differences of max- 
imum likelihood for a given saddle strength follow the 
curve Svi = ± r\/—A for A < 0, independent of which 
velocity difference is being considered (the dashed lines 
for A < in Figs, ^b and ||). Note that the assumption 
that the dominant topology takes the form of a saddle is 
not valid as A nears zero. In this regime, Fig. ^ shows 
that the expected values of the squared strain rate and 
squared vorticity become of comparable magnitude. The 
form of the P{Svi{r)\A) is, therefore, distorted as the 
topological features themselves are stretched (recall Fig. 

Although the distributions of the longitudinal and 
transverse velocity differences are identical across a sad- 
dle, the distribution of these two quantities about a vor- 
tex center differs dramatically. Across an axisymmetric 
vortex center the longitudinal velocity difference is pre- 
cisely zero (since the radial component of the velocity 
itself is zero). For positive values of the Jacobian deter- 
minant, points in the distribution P(A, SvLir)) which lie 
away from 5wl = do so either as the result of first order 
contributions from the strain rate or through higher or- 
der corrections to the flow field. Since the magnitude of 
the longitudinal velocity difference across a vortex center 
is, thus, constrained to be small, it is saddle-like regions 
of the flow which play the dominant role in the wings 
of the longitudinal velocity difference PDFs and, cor- 
respondingly, in the higher order longitudinal structure 
functions. A clear picture of intermittency in 2D must, 
therefore, include an understanding of the formation of 



unusually strong saddle points in addition to coherent 
vortex centers. 

On the other hand, vortex centers do support signifi- 
cant transverse velocity differences. In fact, vortex cen- 
ters with Svt = cannot exist as a consequence of Stoke's 
law: f^vgrfidO = J^tOzd^, where C traces the perime- 
ter of the disk, il, and for A > 0, the sign of vg is the 
same about the entire perimeter. As the magnitude of 
the transverse velocity goes to zero, so must the mag- 
nitude of the vorticity. The shape of the distribution 
is obtained under the assumption that for locations in 
the flow with A 3> the topology is dominated by the 
vorticity. Setting A = (1/4) w^, vg = LurQ, 6vt = 2vg, 
and letting rn = r/2, we find that the transverse veloc- 
ity difference varies as a function of the local Jacobian 
determinant according to Sv^ — zbrVA for A > (the 
dashed line in Fig. ^). The finite spread in the distri- 
bution results both from higher order corrections to the 
shape of the vortex centers as well as from the fact that 
in actuality there are finite contributions from first-order 
saddles, that is, (cr^)/(es) is of order 0.2 for A > 0. 

Because of the larger propensity for strong vortex cen- 
ters to form relative to saddle points of comparable mag- 
nitude (recall Fig. |^a), the maximum transverse veloc- 
ity difference found about centers in the flow is signifi- 
cantly greater than that found about saddles (compare 
Fig. H for A < and A > 0). The wings of the trans- 
verse velocity difference PDF in the enstrophy cascade 
are, therefore, dominated by contributions from vortex 
centers. Because of this clear segregation between the 
structures which play the dominant role in the trans- 
verse and longitudinal velocity differences, it would not 
be surprising if the higher order moments of these two 
quantities differed. 

These rudimentary observations have certainly not ex- 
pended the usefulness of this technique. It is, first of all, 
interesting that the ratio (ct^) /{u!"^) tends toward zero for 
A ^ I. This implies that the strongest vortex centers are 
nearly axisymmetric. Fig. |I|d rather than |l|d. This was 
not seen in the earlier continuously forced experiment 
of Rivera et al. [|T3| . These axisymmetric centers are the 
likely predecessors of the coherent structures found in the 
latter stages of 2D turbulent decay 1 16, 2^. They are sta- 
ble ^ , steady-state solutions of the Euler equa- 
tion in which there is no nonlinear transfer of enstrophy, 
and hence, no enstrophy cascade (27). 

Furthermore, a key feature is still missing from the 
first order distribution, P{6v-!^{r)\A). It is a well know 
result in turbulence theory that the rate of enstrophy 
transfer depends on an odd moment of the longitudinal 



velocity difference (r) — (1/8) yyr'', and, hence, 

on an asymmetry in the distribution of these velocity dif- 
ferences. This asymmetry is missing in Eq. (^, where 
P((5ul|A) = P(— 5ul|A). In fact, any asymmetry is ex- 
plicitly prohibited by incompressibility in the linear field 
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FIG. 4: (a) The expected values of (oj^), filled squares, and 
(cr^), open circles. For positive and negative A different 
topologies tend to be significantly dominant, as is indicated 
by (b), the ratios {uj'^) / {a^) (dashed line) and {a'^) / {uP) (solid 
Une) for negative and positive A respectively. 



approximation. Second order terms, as well, of the form 
(dxxVy) x'^ and {dyyVx)y^ are even under the operations 
X —X and y — > —y so that the contribution to the 
velocity on opposite sides of the structure is of the same 
magnitude and sign. Incompressibility inextricably links 
the remaining dxy terms. The result is that no contribu- 
tion to the longitudinal asymmetry is possible at order 
two. It is necessary to go to at least third order to ex- 
plain the asymmetry in P(Svi^), where the dxxx^x and 
dyyyVy tcrms — aud through incompressibility the dxxyVy 
and d. 



xyy 



terms — result in an asymmetric saddle; that 
is, one in which the magnitude of the velocity in the 
incoming and out-going jets differs. The longitudinal ve- 
locity difference about such a structure is asymmetric, 
whereas, positive and negative transverse velocity incre- 
ments continue to exist with equal likelihood. A system- 
atic correlation between these higher order topological 
structures and the velocity differences is beyond the lim- 
its of the current data set, but future work will explore 



these connections. 
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